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Abstract 

Background: We have used a sensitized ENU mutagenesis screen to produce mouse lines that carry mutations in 
genes required for epigenetic regulation. We call these lines Modifiers of murine metastable epialleles [Mommes). 

Results: We report a basic molecular and phenotypic characterization for twenty of the Momme mouse lines, and 
in each case we also identify the causative mutation. Three of the lines carry a mutation in a novel epigenetic 
modifier, Rearranged L-myc fusion (Rlf), and one gene. Rap-interacting factor 1 (Rifl), has not previously been 
reported to be involved in transcriptional regulation in mammals. Many of the other lines are novel alleles of 
known epigenetic regulators. For two genes, Rlf and Widely-interspaced zinc finger (l/l//z), we describe the first mouse 
mutants. All of the Momme mutants show some degree of homozygous embryonic lethality, emphasizing the 
importance of epigenetic processes. The penetrance of lethality is incomplete in a number of cases. Similarly, 
abnormalities in phenotype seen in the heterozygous individuals of some lines occur with incomplete penetrance. 

Conclusions: Recent advances in sequencing enhance the power of sensitized mutagenesis screens to identify the 
function of previously uncharacterized factors and to discover additional functions for previously characterized 
proteins. The observation of incomplete penetrance of phenotypes in these inbred mutant mice, at various stages 
of development, is of interest. Overall, the Momme collection of mouse mutants provides a valuable resource for 
researchers across many disciplines. 



Background 

Mutagenesis screens for modifiers of position effect 
variegation in Drosophila have played a defining role in 
the development of the field of epigenetics [1]. The 
screens used a fly strain that showed variegated expres- 
sion of the white {w) locus, resulting in red and white 
patches in the eye, as a result of the stochastic establish- 
ment of epigenetic state. The genes identified by these 
screens turn out to have pivotal roles in gene silencing 
[2-4]. We have designed a similar screen in the mouse. 
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using a green fluorescent protein (GFP) transgene 
that shows variegated expression in red blood cells. Off- 
spring of N-ethyl-N-nitrosourea (ENU) -treated males are 
screened for changes in the percentage of erythrocytes 
expressing GFP; this is a screen for dominant effects. 
Each mutant line is named a Modifier of murine metasta- 
ble epiallele Dominant, MommeD [5]. The underlying 
mutations have been identified and published for nine 
lines and the mutations occur in DNA methyltransferases 
{Dnmtl and DnmtSb), chromatin remodelers {SmarcaS 
and Bazlb), a histone deacetylase {Hdacl), a transcrip- 
tional co-repressor {Trim28), a eukaryotic translation ini- 
tiation factor {elFSh) [6-10] and a previously unknown 
gene, Smchdl, now shown to be required for X-inactivation 
in the mouse [11,12]. Recently, Smchdl has been 
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shown to act as a tumor suppressor [13] and mutations 
in SMCHDl have been shown to be tightly associated 
with the human disease facioscapulohumeral dystrophy 
type2 (FSHD2) [14]. 

Results and discussion 

Identification of Momme mutants 

We have now screened a total of approximately 5,000 
Gl offspring, recovered 42 MommeD lines and identified 
the underlying mutations in 29 cases. These lines carry 
mutations in 18 unique genes (Table 1). A number of 
the genes have been hit more than once. Based on their 
effect on GFP expression, 14 lines were classified as sup- 
pressors of variegation, that is, the mutation increased 
the percentage of erythrocytes expressing GFP, and 15 
lines were classified as enhancers of variegation, that is, 
the mutation decreased the percentage of erythrocytes 
expressing GFP (Table 1). Here we report, for the first 
time, the underlying mutations in 20 MommeD lines, 
which carry mutations in 10 unique genes. 

The experimental pipeline for the screen is shown in 
Figure 1. In brief, ENU mutagenesis was carried out as 
described previously [5] and mapping was carried out 
following a G2 backcross to LineSC (a C57BL/6J strain 
carrying the same GFP transgene array at the same loca- 
tion) using traditional microsatellite and SNP genotyping 
or an lUumina GoldenGate SNP genotyping assay. Map- 
ping intervals for the latter group are shown as Manhattan 
plots in Additional file 1. 

To identif)^ the underlying mutations, several different 
strategies have been used. In ten cases, the mutation was 
identified following whole exome deep sequencing of 
DNA from one mutant mouse from the MommeD in- 
bred colony. For these, variants were called within the 
intervals linked to the causative mutation. For some 
lines, variant calling was extended to exome-wide, to as- 
sess for passenger mutations but only a small number 
was found (<10). These did not reflect the typical ENU 
sequence bias (T/A) [30] and probably represent normal 
background mutations. 

In two instances, we have used a custom capture array 
to specifically enrich for DNA from a 4.2 Mbp interval 
previously identified by linkage analysis [7] followed by 
deep sequencing to identify the underlying mutations. In 
the remaining eight cases, Sanger sequencing of the 
exons (including splice sites) of candidate genes in the 
linked intervals was carried out. Details on identification 
of the mutations are provided in Materials and methods. 
In all cases, putative mutations were verified by PGR 
and Sanger sequencing in larger cohorts (at least 100 
mice per mutant line). 

In the two cases in which a custom capture array was 
used, we had an opportunity to obtain an estimate of the 
ENU mutation rate. In one case four mutations were 



found in the 4.2 Mbp captured interval and in the other 
case six mutations were found. Based on these results 
we estimated the ENU mutation rate at approximately 1 
per Mbp. This is consistent with reports by others [31]. 

The heritability of the mutations was tested over at 
least six generations. Flow cytometric expression profiles 
and the percentage of erythrocytes expressing the trans - 
gene for each MommeD (presented here for the first 
time) are shown in Additional file 2. 

MommeDSO mice are haploinsufficient for Widely 
interspaced zinc finger (Wiz) 

MommeDSO was classified as an enhancer of variegation. 
We mapped the MommeDSO mutation to a 1.9 Mbp 
interval on chromosome 17 encompassing 48 genes 
(Additional file 3). Exome deep sequencing identified a 
single base deletion in exon 5 of the Wiz gene, causing a 
frame shift mutation that is predicted to introduce a pre- 
mature stop codon (Figure 2a and Table 1). No other pu- 
tative ENU variants were identified in the linked region. 

Wiz contains six Kruppel-type zinc finger motifs in a 
widely interspaced manner and has been reported to associ- 
ate with the histone H3 lysine 9 (H3K9) methyltransferases 
G9a and GLP in cell lines [32,33]. No mouse mutants for 
Wiz have so far been described. Western blot analysis using 
an anti-Wiz antibody detected two bands (about 120 and 
130 kDa), as expected, in wild-type embryos (Figure 2a). 
Embryos heterozygous for the MommeDSO mutation 
showed approximately half the amount of Wiz protein and 
no Wiz protein was detected in homozygotes (Figure 2a). 
We conclude that MommeDSO mice are haploinsufficient 
for Wiz. We refer to this allele as ^7/^^^^^^^^. 

MommeDIS carries a mutation in Rap interacting factor 1 
(Rim 

MommeDIS was classified as a suppressor of variegation 
and candidate gene sequencing revealed that MommeDIS 
mice carry a nonsense mutation in exon 29 of the Rifl 
gene (Figure 2b and Table 1). Two other candidate genes 
{Epc2 and MbdS) were sequenced and no mutations were 
identified in these genes. Western blot analysis showed 
approximately half the amount of full-length Rifl protein 
in MommeDIS heterozygotes and no full-length Rifl 
protein in a homozygote (Figure 2b). Furthermore, a Rifl 
gene trap allele (Materials and methods) had a similar ef- 
fect on transgene expression as that observed with the 
MommeDIS mutation, increasing the percentage of ex- 
pressing cells in mice heterozygous for the gene-trap allele 
(Additional file 4). 

We conclude that the MommeDIS allele is the result 
of a mutation in Rifl that we refer to as Rifi^ommeDis ^ 
Rifl has recently been shown to be involved in DNA 
double strand break repair [34-37] and replication 



Table 1 MommeD mutants, causative mutations and disease association 



Name 


Effect on variegation 


Gene 


Nature of mutation 


Homozygous phenotype 


References 


Human homolog 


Disease association 


MommeD 1 


Suppressor 


Smchd! 


C->T exon 23; introduces Stop 


Null 


[5,11] 


SMCHDl 


FSHD2 [14] 


MommeD2 


Suppressor 


Dnmt! 


C->A exon 25; T812K 


Null 


[6] 


DNMTl 


Schizophrenia, breast and prostate 
cancer l 1 5- 1 yj 


MommeD4 


Er^hancer 


5m areas 


T->A exon 12;W520R 


Hypomorphic? 


[5,6] 


SMARCA5 


Acute myeloid leukemia [20] 


MommeDS 


Enhancer 


Hdac! 


7 bp deletion in exon 13; Frameshift 


Null 


[5,7] 


HDACl 


Schizophrenia, neural development [2] 


MommeDS 


Enhancer 


RIf 


G->T exon 8;C1558F 


Hypomorphic 


This study 


RLE 




MommeD9 


Enhancer 


Trim28 


T->C splice donor site of intron 13 


Null 


[8] 


TRIM28 




MommeD! 0 


Enhancer 


Baz!b 


T->G exon 7; L733R 


Hypomorphic 


[7] 


BAZIB 


Williams-Beuren syndrome [21] 


MommeD! 2 


Enhancer 


elF3h 


T->A - 10 bp before splice acceptor 
site of intron 4 


Null 


[9] 


elF3H 




MommeD! 3 


Suppressor 


beiuui 


A->G exon 20; results in splicing defect 


Ml ill 
NUN 


This study 


OL 1 Ud I 


Melanoma [22] 


IVIOllUlieU 1 H 


Suppressor 


UtItllLju 


T->C splice acceptor site of intron 12 


Hypomorphic 


L 1 uj 


UINIVI 1 JD 


ICE syndrome [23] 


MommeD! 6 


Enhancer 


Dazi 0 


1 ->L. exon z, L/jr 


Hypomorphic? 


This study 


RA71 R 
dAA I D 


Williams-Beuren syndrome [21] 


MommeD!? 


Suppressor 


Setdb! 


1 ->L exon z 1 , V 1 z4oA 


Hypomorphic 


This study 


SETDBl 


Melanoma [22] 


MommeD! 8 


Suppressor 


Rif! 


C->T exon 29;Q1669 Stop 


Null 


This study 


RIFl 


Breast cancer [24] 


MommeD! 9 


Suppressor 


5m a red 


T->G splice acceptor site of intron 10 


Null 


This study 


SMARCCl 


Colorectal cancer [25,26] 


MommeD23 


Suppressor 


Smchd! 


A->T exon 1 2; R498 Stop 


Null? 


This study 


SMCHDl 


FSHD2 [14] 


MommeD27 


Suppressor 


Pbrm! 


A->G exon 17; Y733C 


Hypomorphic? 


This study 


PBRMl 


Renal cancer [27] 


MommeD28 


Enhancer 


RIf 


A->G splice acceptor site of intron 4 


Null 


This study 


RLE 




MommeD30 


Enhancer 


Wiz 


1 bp deletion in exon 5; Frameshift at 
amino acid 553 


Null 


This study 


WIZ 




MommeD3 ! 


Enhancer 


Trim 28 


T->A exon 3;C178S 


Null? 


This study 


TRIM28 




MommeD32 


Suppressor 


Dnmt! 


T->C exon 29; L1045P 


Null? 


This study 


DNMTl 


Schizophrenia, breast and prostate 
cancer l i d- i yj 


MommeD33 


Suppressor 


buvjyn 1 


A->G splice donor site of intron 1 


Null 


This study 


CI IX/QOU1 1 


Lupus, retinoblastoma [2] 


MommeD34 


Enhancer 


DIf 
KIT 


L->A exon /, Lojj otop 


Ml ill 
IMUII 


This study 


Dl C 




MommeD35 


Enhancer 


SmarcaS 


A-> (j exon y, i\io4ib 


Hypomorphic? 


This study 


A or" A c 
oIVIAKLAd 


Acute myeloid leukemia [20] 


MommeD36 


Suppressor 


Smchd! 


C->T exon 42; Q1732 Stop 


Null? 


This study 


SMCHDl 


FSHD2 [14] 


MommeD37 


Enhancer 


SmarcaS 


T->C exon 13; L565P 


Null? 


This study 


SMARCA5 


Acute myeloid leukemia [20] 


MommeD38 


Enhancer 


elF3h 


G->A exon 7; R291 Stop 


Null 


[9] 


elF3H 




MommeD39 


Suppressor 


Smarca4 


G->A splice donor site of intron 20 


Null 


This study 


SMARCA4 


Coffin-Siris syndrome [28] 


MommeD40 


Suppressor 


Uhrf! 


T->A exon 17; Y778 Stop 


Null 


This study 


UHRFl 




MommeD42 


Enhancer 


Brd! 


T->A exon 11; C411 Stop 


Null 


This study 


BRDl 


Schizophrenia, bipolar affective 
disorder [29] 



5 i. 

g a. 

O rt> 

^§ 

O C3D 

§! 

.S" 
O ^ 

^ NJ 
O 

^ 
VO - 

VO 



In some cases the mutation can be defined as hypomorphic or null and in other cases this designation is less certain and we have added a question mark. 
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A male FVB/N J with 
a variegated and epi- 
genetically sensitive 
GFP transgene is 
injected with the 
mutagen END. 



Backcross to 
C57BL/6 for 
two generat- 
ions allows 
for linkage 
analysis. 




The offspring are 
screened for changes : 
in GFP expression j 4x 
using a drop of blood. ; 



C57BL/6 r 



60 
50 
C 40 
D 30 
O20 
^ 10 



~ Mutant 
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. GFP* 



1 



Log fluorescence 

Mice carrying mutations 
that shift the percentage 
of GFP expressing 
erythrocytes are tracked 
by FACS. 



Iff 



52 ^ 



Backcross 



. /VS .. 



lllumina GoldenGate Genotyping 
Assay (12 mutants, 12 wild types), 
and fine-mapping, link the causative 
mutation to a genomic interval. 



ENU-induced point mutations are identified 
using whole exome sequencing, or gene 
candidate based sequencing. 



Figure 1 An ENU mutagenesis/whole exome deep sequencing pipeline enables rapid identification of causative mutations in mice with 
defects in epigenetic gene silencing. A schematic overview of tlie ENU mutagenesis gene discovery pipeline is presented. Tine major 
components of tine screen are described in tine figure. Briefly, male FVB/NJ mice carrying an epigenetically sensitive GFP transgene {LineS) were 
treated with ENU and mated with female LineS mice. Offspring were screened for a shift in the percentage of GFP expressing erythrocytes using 
flow cytometry. Putative mutants were mated with LineS mice for four generations to test for heritability and reproducibility of the GFP 
expression, and to reduce the number of non-causative ENU mutations within the genomes. Linkage analysis was carried out on the offspring 
from two generations of backcrosses between putative mutants and C57BL/6J mice, which also carry the GFP transgene {LineSQ. Causative 
mutations were identified by whole exome deep sequencing, or gene candidate sequencing on individuals that had been maintained on the 
Lines background for at least seven generations. FACS, fluorescence-activated cell sorting. 



timing in mammals [38,39], both processes involving sig- 
nificant chromatin reorganization. Our findings suggest 
a role for Rifl in transcription. 

MommeDS, MommeD28 and MommeD34 carry mutations 
in Rearranged L-myc fusion (RIf), a gene about which little 
is known 

Three enhancers of variegation, MommeDS, MommeD28 
and MommeD34y were found to carry mutations in Rlf, 
Whole interval capture revealed that MommeDS and 
MommeD34 mice have mutations in exon 8 and exon 7, 
respectively. MommeDS is a missense mutation and 
MommeD34 is a nonsense mutation (Figure 2c and 
Table 1). The MommeD2S mutation was identified by can- 
didate sequencing and affects the splice acceptor site of in- 
tron 4 (Figure 2c and Table 1). Western blot analysis using 
two independent anti-Rlf antibodies (one polyclonal and 
one monoclonal) detected an approximately 280 kDa 



band in protein lysates made from wild-type embryos 
(Figure 2d). Although the band detected on the western 
blots is substantially higher than the predicted molecular 
weight of 218 kDa (NCBI Mouse Build 37), this could be 
due to post-translational modifications to the Rlf protein. 
We concluded that it is likely that this band is Rlf. 
Reduced amounts (MommeDS) or no {MommeD2S and 
MommeD34) Rlf protein was detected in homozygous em- 
bryos (Figure 2d). Together, these results suggest that 
MommeDS is a hypomorphic allele and that MommeD2S 
and MommeD34 are null alleles. 

Rlf encodes a protein predicted to contain 16 widely 
spaced zinc finger domains, which has led to the sugges- 
tion that it has a role in transcriptional regulation [40]. 
Rlf is conserved across vertebrates. To determine the 
localization of Rlf in the cell, we carried out cell frac- 
tionation experiments. Rlf protein was specifically 
detected in the nuclear fraction of HeLa cells (Figure 2e). 



Daxinger et al. Genome Biology 201 3, 14:R96 Page 5 of 1 7 

http://genonnebiology.conn/201 3/1 4/9/R96 



MommeD30 

F ram ash if t at aa 553 



= Zn Finger Domain 



MommeD30 AGRPGNQELGRPKF 
Mouse AGRPGKPGAGPAQV 
Human AGRPGKPGAGPTQV 



Wiz 



956 aa 



MommeD18 

Q 1669 X 



Rif1_N 



Rif1 



2426 aa 



MouseRLLYKQDNDSN 
Human RL - - K R D I DSN 



+/+ 



Wiz Protein 

D30/+ 



Wiz 



Y-tubulin 



D30/D30 (kDa) 
150 



+/+ 



Rifl Protein 
D18/+ 



100 
t 50 



MommeD28 

Splicing Defect 



MommeD34 

C355X 



■■ ■■■■■ 



RIf 



= Zn Finger Domain 

■ ■■ 



MommeDS 

C 1558 F 

I 



3IED: 



Mouse TTTCTAG 
Human TTTCTAG 
Opossum T C C I TAG 
Stickleback 



\ 

VN KL I 
VN KL I 

A-v N R r L 



Intron 4 Exon 5 



r 



QQHLFCL I RV I 
QQHLFCLI RV I 
QQHLFCLI RVV 
V HHLLCLVRV I 



Rif Protein 



( 5^ \ 

Mouse H TQYPCM VQGC 
Human H TQYPCM VQGC 
Opossum H TQYPCM QGC 
Stickleback IRAAYPCMVQDC 



D18/D18 (kDa) 




1918 aa 



+/+ D8/D8 

RIf (polyclonal) 
RIf (monoclonal) — 
T^tubulin 



+/+ D28/D28 




+/+ D34/D34 (kPa) 

- 250 



N (kDa) 




RIf 
GAPDH 
SMARCA5 



250 
37 
150 



f 




+/+ 

n=94 

D8/+ 

n=95 








=5 




:::::::::::: 82 :::::::::::: 





9 



12 




63 




1 


—I 1 1 1 



0 20 40 60 80 

Mouse Coat Colour 

(% of cohort) 

**Chi-test< 0.0001 

Figure 2 (See legend on next page.) 



100 



Daxinger et al. Genome Biology 2013, 14:R96 
http://genonnebiology.conn/201 3/1 4/9/R96 



Page 6 of 1 7 



(See figure on previous page.) 

Figure 2 Causative mutations in MommeDSO, MommeDIS, MommeDS, IVIommeD28 and IVIommeD34. (a) MommeDSO carries a 1 bp 
deletion in l/l//ztliat leads to a frame-shift at amino acid 553. Western blot analysis of embryo heads at 12.5 days post-coitum (dpc) shows 
reduced levels of Wiz protein in \/\/iz^^^^^^^o heterozygotes and no Wiz protein in homozygotes. Wiz protein was detected at approximately 120 
and approximately 130 kDa. Each track represents a different animal. Tubulin was used as a loading control, (b) MommeDIS harbors a point 
mutation in Rifl that introduces a premature stop codon at amino acid 1,669. Western analysis of testes tissue shows that f^jf^ommeDis 
heterozygotes have a reduced dosage of Rifl. No Rifl protein was detected in a homozygote. Each track represents a different animal. Rifl 
protein was detected at approximately 260 kDa and Tubulin was used as a loading control, (c) ENU point mutations in MommeDS, MommeD2S 
and MommeD34 occur in conserved regions of the RIf protein. Asterisks indicate mutation, (d) Western blotting of RIf in embryo head lysates 
from 14.5 dpc Rlf^^, RlfMomme^^^^^, Rlf^^^^^^ and Rif^^^^^^ revealed greatly reduced amounts of RIf protein. Each track represents a different 
animal. RIf protein was detected at approximately 280 kDa. Tubulin was used as a loading control, (e) Total, cytoplasmic and nuclear fractions of 
HeLa cells were isolated and protein concentration quantified. Equal amounts of each fraction were immunoblotted with anti-RIf, GAPDH 
(cytoplasmic marker) and SMARCA5 (nuclear marker) antibodies, revealing nuclear localization of RLE. RIf protein was detected at approximately 
280 kDa, GAPDH protein at 37 kDa and SMARCA5 at approximately 120 kDa. (f) Coat color of offspring carrying the allele produced from a 
^ifAommeDs heterozygous female crossed to a pseudoagouti /\^^ male. Rif^"^^"^^^ heterozygous offspring showed a shift in coat color towards 
pseudoagouti compared to wild-type littermates. 



To our knowledge, we have identified the first mouse 
mutants for RIf and hereafter refer to these alleles as 

We have identified RIf from a screen that is based on 
the expression of a multi-copy variegating transgene. To 
test if RIf is required for the establishment of epigenetic 
state at an endogenous locus, we studied the effect of 
j^l^ommeDs gxpression of agouti viable yellow {A'^), 
Ai^ is a single copy gene that, like the transgene, is 
known to be sensitive to epigenetic state [41]. The A^^ 
allele is the result of an intracisternal-A-particle (lAP) 
retrotransposon insertion upstream of the agouti gene 
and isogenic mice carrying the A^^ allele can be yellow, 
mottled or pseudoagouti (dark brown). Coat color has 
been shown to correlate with the level of DNA methyla- 
tion at the long- terminal repeat (LTR) promoter [42]. It 
has been shown that haploinsufficiency for modifiers of 
epigenetic reprogramming can alter the ratio of the dif- 
ferent coat colors [5,6,43,44]. We crossed females het- 
erozygous for j^i^o^^^^^ ^ith pseudoagouti A^^/a mice 
and scored the offspring for coat color. Offspring that 
inherited the Rlf^ommeDs ^^^^^ ^^^^ ^^^^ ^.^^^^^ 

pseudoagouti than their wild-type littermates, that is, 
haploinsufficiency for RIf increased the probability of si- 
lencing at the A^^ locus (Figure 2f). This is consistent 
with RIf being an enhancer of variegation and suggests 
that RIf hsiS a general role in epigenetic regulation. 

Novel alleles of known epigenetic regulators 

In addition to the genes described above, we have pro- 
duced 15 MommeD lines that carry mutations in genes 
already known to be involved in epigenetic regulation in 
the mouse. Exome deep sequencing of three suppressors of 
variegation, MommeD 19, MommeD27 and MommeD39, 
identified mutations in the chromatin remodelers SWI/ 
SNF related, matrix associated, actin dependent regulator 
of chromatin, subfamily c, member 1 {Smarccl), Polybromo 
1 {Pbrml) and SWI/SNF related, matrix associated, actin 
dependent regulator of chromatin, subfamily a, member 4 



MommeD 1 9 



carries 



{Smarca4), respectively (Table 1). Smarccl 
a mutation in the splice acceptor site of intron 10 predicted 
to lead to an in-frame premature stop codon and likely 
nonsense-mediated mRNA decay. Quantitative real-time 
RT-PCR and western blot analysis showed that Smarccl 
mRNA and Smarccl protein levels were reduced in hetero- 
zygous animals (Figure 3a). The Pbrml^'''^'^'^^^ mutation 
is a missense mutation in exon 17. Western blot analysis 
revealed reduced levels of Pbrml protein in homozygotes 
(Figure 3b). Smarca4^'''^'^'^^^ has a mutation in the splice 
donor site of intron 20. cDNA sequencing revealed an ex- 
tended exonic sequence (data not shown) that is predicted 
to lead to an in-frame premature stop codon and 
nonsense-mediated mRNA decay. Quantitative real-time 
RT-PCR analysis showed reduced Smarca4 mRNA in het- 
erozygotes (Figure 3c). 

MommeD42 was identified as an enhancer of variega- 
tion and carries a nonsense mutation in exon 11 of 
Bromodomain containing 1 {Brdl; Figure 3d and Table 1). 
Brdl has recently been reported to form a complex with 
the histone acetyltransferase HBOl and is required for 
the transcriptional activation of the erythroid-specific 
regulator genes in fetal liver [45]. 

Three suppressors of variegation were found to carry 
mutations in H3K9 methyltransferases. MommeDIS and 
MommeDl? carry mutations in Setdbl and MommeDSS 
has a mutation in Suv39hl (Table 1). The Setdbl'^''^'^''''^ 
mutation introduces a novel splice donor site leading to 
an in-frame premature stop codon in exon 20 (Figure 3e). 
cDNA analysis confirmed the presence of the predicted 
transcript (data not shown). Setdbl^'''^'^^^^^ is a missense 
mutation in exon 21 of Setdbl (Figure 3e). RNA analysis 
indicates reduced levels of Setdbl mRNA in the 
MommeDIS mutants whereas Setdbl mRNA levels were 
not significantly different in the MommeDl? mutants 
(Figure 3e). The Suv39hl^''^'^'''^^ mutation is at the 
splice donor side of intron 20. This mutation is predicted 
to lead to a premature stop codon. Consistent with this, 
western blot analysis showed that the level of Suv39hl 
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(See figure on previous page.) 

Figure 3 Causative mutations in MommeD19, MommeD27, MommeD39, MommeD42, MommeDIS, MommeDSS and IVIommeD40. 

(a) MommeD19 carries a mutation at a splice site of Smorccl. Real-time RT-PCR and immunoblotting of embryos at 12.5 days post-coitum (dpc) 
showed reduced levels of Smarccl mRNA and Smarccl protein in heterozygotes (n > 5 mice). Each track represents a different animal. Smarccl 
protein was detected at approximately 160 kDa. (b) MommeD27 harbors a mutation in a bromodomain of Pbrml. Western blot analysis of 14.5 
dpc embryo heads showed reduced levels of Pbrml in homozygotes. Each track represents a different animal. Pbrml was detected at 
approximately 190 kDa. (c) MommeD39 carries a mutation at a splice site of 5marca4. Real-time RT-PCR showed reduced Smarca4 mRNA in testes 
of heterozygotes (n > 4 mice). cDNA sequencing revealed that the mutation in Smorco4'^°^^^^^^ results in use of an alternative splice donor site, 
(d) The mutation in MommeD42 introduces a premature stop codon at amino acid 41 1 of the Brdl protein, (e) Mutations in MommeDIS and 
MommeDU occur in the conserved SET domain of Setdbl. Real-time RT-PCR of Setdbl mRNA from testes of heterozygotes and age-matched 
wild types (n = 4 mice). cDNA analysis revealed that the Setdbl'^°^^^^^^ allele is associated with the use of an alternative splice donor site in 
exon 20, leading to a 62 bp truncation, (f) MommeDSS carries a mutation at a splice site in SuvS9hl. Western blot analysis of Suv39hl in adult 
thymus showed reduced Suv39hl in hemizygous mutant males. Each track represents a different animal. Suv39hl protein was detected at 48 
kDa. (g) The mutation in MommeD40 introduces a premature stop codon at amino acid 778 of Uhrfl. Western blot analysis of Uhrfl revealed 
greatly reduced levels in 9.5 dpc embryos homozygous for the iii-ii^f]^ommeD4o i^y^g^jon. Each track represents a different animal. Uhrfl protein 
was detected at approximately 90 kDa. Error bars indicate ± standard error of the mean. N.s., not significant. Asterisks indicate mutation. 



protein was greatly reduced in hemizygous mutant males 
(Suv39hl is on the X chromosome; Figure 3f). 

MommeD40, a suppressor of variegation, carries a non- 
sense mutation in exon 17 of the Uhrfl gene (Figure 3g 
and Table 1). Western blot analysis revealed greatly re- 
duced levels of Uhrfl protein in homozygous embryos 
(Figure 3g). Uhrfl is known to associate with the mainten- 
ance DNA methyltransferase Dnmtl [46] and has recently 
been shown to link DNA methylation and H3K9 methyla- 
tion in human cells [47]. 

Seven MommeD lines are new mutant alleles of 
genes already identified in the screen and we have des- 
ignated these lines Bazlb'^''^'^''''^ Smchdl^''^'^'''^^ 
Trim28'^'''^'^'''^\ Dnmtl'^'''^'^'''^^ SmarcaS'^'''^'^'''^' , 
Smchdl^'''^'^'^^^ and SmarcaS'^'''^'^'^^^ (Table 1; 
Additional file 5). 

Approximately 60 Suppressor of variegation {Su(var)) 
and 25 Enhancer of variegation {E(var)) genes were iden- 
tified in the Drosophila screens [2]. We have now docu- 
mented 12 Su(var) and 8 E(var) genes in the mouse and 
these encode some of the key factors for DNA methyla- 
tion and H3K9 methylation, as well as a large group of 
chromatin remodelers (Figure 4). Some of these genes, 
such as Dnmtl y Suv39h and Trim28, have been identi- 
fied as core heterochromatin factors, whereas others, 
such as SmarcaS and Hdacl, have been associated with 
gene silencing at euchromatic positions [2]. This is con- 
sistent with the variegating phenotype of the transgene. 
Interestingly, we have not yet recovered any members of 
the RNA interference/piwi-interacting RNA (piRNA) 
pathway or Polycomb and trithorax groups of proteins 
but the screen has not reached saturation. 

Momme genes are required for normal embryonic 
development 

In most of the MommeD lines, heterozygous mice were 
observed at expected ratios at weaning and did not show 
any obvious phenotypes. To determine the viability of 
homozygotes, we performed embryonic dissections 



following intercrosses. We have focused on those lines 
carrying mutations in genes about which there is no or 
few data available in the literature. 

Following \)^i^ommeD30 i^tercrosses, no homozygous 
animals were recovered at weaning out of 97 progeny 
(Figure 5a). To determine when w^/^^^^^^^^ homozy- 
gotes died, embryos from intercross matings were 
obtained at different stages of development and geno- 
typed. Viable homozygotes were recovered at the 
expected ratios at 10.5 days post-coitum (dpc). However, 
these embryos appeared smaller than their wild-type and 
heterozygous littermates. At 12.5 dpc homozygotes were 
present in numbers less than expected (15% as opposed 
to 25%) and at 14.5 dpc no viable homozygous embryos 
were obtained (Figure 5a). We conclude from these data 
that embryonic death in mice homozygous for 
^^^ommeDso ^^^^^^ between 10.5 and 12.5 dpc. 

Rifl^ommeDis i^tercrosses produced wild-type and het- 
erozygous offspring at the expected ratios and occasional 
homozygous individuals (3 out of 155 offspring) at 
weaning (Figure 5b). All surviving homozygotes were 
males and two, out of three tested, were infertile. Timed 
matings revealed that the majority of the homozygotes 
died during the second half of gestation (Figure 5b). This 
is consistent with the literature [48]. 

i?//' intercrosses produced homozygous individuals that 
survived to weaning in less than expected numbers for 

Rlj^ommeDS^ j^ljMommeD28 j^yMommeD34 (pig^^.^ 5^)^ 

j^l^ommeD34 homozygotes died in the first week after 
birth, for unknown reasons (Figure 5c). A similar finding 
was made for Rlj^ommeD28 homozygotes; most died in 
the first week after birth (data not shown). However, in 
both lines some homozygotes did survive to adulthood 
and were fertile (data not shown). Embryonic dissections 

indicated that, for the RlfommeD28 j^i^ommeD34 

lines, late gestation homozygous embryos weighed sig- 
nificantly less than their heterozygous or wild-type litter- 
mates (T-test, P < 0.05; Figure 5c). RlfommeDs b^^aves 
like a hypomorphic allele, with approximately half the 
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to be involved in transgene silencing in the mouse. Factors are grouped by their broad mechanistic role in epigenetic regulation. TF, transcription 
factor; Me, methyl. Figure adapted from [2]. 



expected number of homozygotes present at weaning 
(Figure 5c). Those surviving were smaller and were fer- 
tile (data not shown). 

One advantage of ENU mutagenesis is that it produces 
point mutations, increasing the likelihood of hypomorphic 
mutations. We have reported hypomorphic alleles for 
SmarcaS'^'''^'^'''^ and DnmtSb'^'''^'^''''^ previously [6,10]. 
In addition to the Rlf^^^'^^ allele described above, 
Setdb'^'''^'^''''^ behaves like a hypomorphic allele for 
the following reasons. Following intercrosses, some 
Setdbl^^^^^^^^ homozygotes were present at weaning 
(Figure 5d) and some survived to adulthood but were 
smaller than their wild-type and heterozygous littermates 
and showed decreased fertility (data not shown). Mice 
homozygous for a null allele of Setdbl have been reported 
to die at the early post-implantation stage [49]. Consistent 
with this, timed matings with Setdbl^''^'^''''^ heterozy- 
gous mice produced no viable homozygous offspring at 
9.5 dpc (Figure 5d). To our knowledge, Setdbl'^'''^'^''''^ is 
the first hypomorphic allele for Setdbl, 



In two cases, Smarcd 



MommeDl 9 



and Smarca4 



\MommeD39 



heterozygous mutants showed signs of stochastic death 
prior to weaning. Following Smarcd 



MommeD 19 



inter- 



crosses, no viable homozygotes were obtained at 
weaning and timed matings revealed that homozygous 
lethality occurred prior to mid-gestation, as expected 
based on the literature [50] (Figure 5e). At weaning the 
percentage of heterozygotes to wild types (47% and 53%, 
respectively) was less than that expected (66% and 33%, 
respectively). Maintenance of the colony, which involved 
heterozygous to wild type crosses, also produced fewer 
mutants than expected at weaning (34% as opposed to 
50%) (Figure 5e). Further investigation suggested that 
this stochastic death of heterozygotes was occurring 
after 17.5 dpc and before 1 week postnatal (Figure 5e). 



At 17.5 dpc, some heterozygotes had exencephaly, 
although this could not account for all the perinatal le- 
thality observed (Figure 5e). The average weight of het- 
erozygous mice at E17.5 and at 1 week postnatal was 
reduced (T-test, P < 0.0001) and greater variance (F-test, 
P < 0.05) was observed among heterozygotes than 
among wild- type littermates (Figure 5e). 

Maintenance of the Smarca4^'"^'^''^^^ colony (hetero- 
zygote to wild type crosses) also revealed non- 
Mendelian ratios; heterozygotes were present in less 
than expected numbers at weaning (18% as opposed to 
50%) and were smaller than their wild-type littermates 
(Figure 5f). Timed matings showed that the ratios were 
close to Mendelian prior to birth (Figure 5f). Body 
weights of heterozygotes were lower than those of wild- 
type littermates at weaning (T-test, P < 0.0001) and, 
again, the variance was greater in the heterozygotes 
(F-test, P < 0.0001) (Figure 5f). Similar stochastic death 
of heterozygotes has been described by others for mice 
with knock-out alleles for Smarcd and Smarca4 [50,51]. 
In these cases, however, the genetic heterogeneity among 
the mice could have been an underlying factor. 

We have previously reported increased variance in 
body weight and behavioral responses in adult mice het- 
erozygous for Trim28 [8] and others have made similar 
findings with respect to Trim28 during the oocyte to 
embryo transition [52]. We have also reported that re- 
duced levels of Dnmt3a are associated with an increase 
in body weight variance in adults [8]. These observa- 
tions, in the context of inbred strains reared in con- 
trolled environments, are consistent with a role for the 
genes identified in this screen in canalization, a term 
coined by Waddington to describe phenotypic robust- 
ness during development. Exactly how important prob- 
abilistic developmental events are in determining 
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(See figure on previous page.) 

Figure 5 MommeD mutants show abnormal embryonic development. Timed matings and intercrosses of Mom meD mutants, 
(a) i/i///^ommeD3o j^j^^. ^^^^ s|-iow tlie number of mice observed (and in bracl<ets tine percentage) at 10.5 days post-coitum (dpc), 12.5 dpc, 14.5 
dpc and at weaning. Representative embryos are sliown. (b) f^jf]^^^^^^^^ m\ce; data sliow tine number of mice observed (and in bracl<ets tine 
percentage) at 10.5 to 12.5 dpc, 14.5 dpc and at weaning, na, not applicable, (c) /?/f mutants; data show the number of mice observed (and in 
brackets the percentage). Embryonic weights were measured from intercrosses of f^if^°^^^^^^ and f^if^^^^^^^^ m\ce. Homozygous embryos from 
both had a significant increase in weight variation at 17.5 dpc or 16.5 dpc, respectively. Weights for each litter were normalized to the average 
weight of wild-type embryos in that litter. Each data point represents an individual. Homozygotes were smaller than wild-type littermates at one 
week after birth, ns, not significant, (d) Setdbl^"'^'^''^^^ and Setdbl^"'^'^''^^^ mice; data show the number of mice observed (and in brackets the 
percentage), (e) Smorccl^"^"^^^^^ and Smorccl^"^^^^^^ mated to wild-type mice revealed incomplete penetrance of heterozygous lethality after 
birth. Data show the number of mice (and in brackets the percentage). Weights of 17.5 dpc Smorccl'^°^^^^^^ embryos and pups were, on 
average, less than that of wild-type littermates (T-test, P < 0.0001) and showed greater variation (F-test, P < 0.05). Weights in each litter were 
normalized to the average weight of wild-type embryos in that litter. Each data point represents an individual, (f) Smorco4^°"^^^^^^ heterozygotes 
showed reduced viability. Data show the number of mice and in brackets the percentage. Heterozygotes were smaller than their wild-type 
littermates. Smarc4^°'^'^^^^^^^ weighed less (T-test, P < 0.0001) and weights were more variable (F-test, P<0.0001) than wild types. Weights in 
each litter were normalized to the average weight of wild-type embryos in that litter. Each data point represents an individual. 



phenotype is not yet clear [53]. We have so far reported 
only three measures; the penetrance of homozygous or 
heterozygous death in the colony, the body weight 
of heterozygotes or homozygotes (at various ages) in 
the colony and behavioral responses. More extensive 
phenotyping of these inbred MommeD lines should 
enable us to gain a better understanding of this 
phenomenon. 

In the remaining MommeD lines examined, 
Bazlb^'''^'^'^^^, Smchdl^'''^'^'^^^ Pbrml^'''^'^'^^^ , 
Trim28^'''^'^'^^\ Dnmtl^'''^'^'^^^, Suv39hl^'''^'^'^^^ 
SmarcaS'^''^'^'''^^ , Smchdl'^''^'^'''^^ SmarcaS'^'''^'^'''^^ , 
Uhrfl'^'^'^'''^^ and Brdl^''^'^'''^' , we observed homo- 
zygous embryonic death that is similar to that reported 
for knock-out alleles of these genes (Additional file 6). 



DNA methylation levels at the transgene are increased in 
Rlf mutants 

Changes in transcription can correlate with changes 
in DNA methylation at transgenes and metastable 
epialleles [42,54-56]. We have shown previously that 
changes in the percentage of red blood cells expressing 
GFP can be accompanied by changes in DNA methyla- 
tion at the transgene [7,10]. Using bisulfite sequencing, 
we investigated DNA methylation at the HS-40 enhan- 
cer region of the transgene in MommeD lines. We used 
adult spleens from wild-type mice and mice heterozy- 
gous for a MommeD mutation. Consistent with our pre- 
vious reports, in LineS wild-type mice around 60% 
of the CpGs in the HS-40 element were methylated 
(Figure 6a). Mice heterozygous for Dnmtl^'''^'^'^^^ sug- 
gested a decrease in CpG methylation (around 50%; 
Figure 6a). This is consistent with an increase in expres- 
sion of the GFP transgene in these mutants and with its 
role as the maintenance DNA methyltransferase. In 
mice heterozygous for the Wiz and Rifl mutations DNA 
methylation patterns at the HS-40 region were un- 
affected (Figure 6b). 



In the case of the Rlf alleles, we detected higher lev- 
els of DNA methylation at the HS-40 element in 
adult spleens in Rlfommeos (320/0, T-test <0.05) and 
RlfommeD34 ^gp/^^ homozygotes Compared 

to wild types (67%) (Figure 6d). A similar trend (not sta- 
tistically significant) was observed in the Rlf^^^^^^ 
and the Rif^^^^^^^^ heterozygotes compared with that 
seen in LineS wild-type (Figure 6a,c). In the case of 
j^ljMommeD28 ^ where adult homozygotes are extremely 
rare, we analyzed DNA methylation patterns in 16.5 dpc 
embryos and made a similar finding. The HS-40 elem- 
ent was more methylated in Rif^^^^^^^^ homozygotes 
(79%, T-test <0.05) than it was in wild-type controls 
(51%) (Figure 6d). This is the first time that we have ob- 
served hypermethylation of the transgene in a MommeD 
line. An increase in DNA methylation at the transgene 
locus in the i?//' mutants is consistent with the decreased 
expression of the transgene in these mice. 



Conclusions 

Our results suggest that transgene silencing in the 
mouse acts through a mechanism common to trans- 
poson silencing, X-inactivation and imprinting. We pro- 
vide evidence that a novel gene, Rlf is involved in this 
process. Depletion of i?//' leads to DNA hypermethylation 
at the transgene. How this is achieved is currently un- 
known and requires further investigation. We have iso- 
lated the first mouse mutants for Rlf and Wiz and 
shown that the genes are required for normal embryonic 
development. Notably, the human homologs of many of 
the genes recovered from our screen have been found to 
be associated with human diseases, in some cases identi- 
fied in family studies and in other cases suggested by 
genome- wide association studies (Table 1). 

Since the mouse is often used as a disease model, we 
anticipate that our collection of Momme mutants 
will provide a valuable resource for researchers across 
many disciplines. 
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(See figure on previous page.) 

Figure 6 Bisulfite sequencing of the HS-40 enhancer region of the GFP transgene. (a) DNA metliylation levels were reduced in 
Dn/T]t/^°"^"^^^"^^ heterozygous mice compared to wild-type mice in adult spleens (n = 4 wild-type mice and n = 2 /3/^/77t/^ommeD32 p^y^gp^ mice). 
Filled circles represent methylated CpG sites, (b) DNA methylation levels in adult spleens were unaffected in mice heterozygous for i/i///^ommeD3o 

and RjfMommeDlS ^ ^ ^■^ommeD30 n = 2 for Rjfl^ommeDlSy methylation levels in adult spleens of RjfommeDS^ ^if^ommeD28 

^l^ommeD34 heterozygotes were increased compared to wild-type mice (n = 2 mice for all genotypes), (d) /?/f homozygotes showed significantly 
increased DNA methylation (J-test, P < 0.05) in adult male spleen (Rlf^^^^^Ds ^jfAommeD34^ embryo [Rif^°^^^D28^ compared to 

wild types (n = 3 for all genotypes). 



Material and methods 

Generation and screening of mutant mice 

Procedures were approved by the Animal Ethics Commit- 
tee of the QIMR Berghofer Medical Research Institute. 
The ENU screen was carried out in the FVB/NJ inbred 
transgenic line, LineS, which is homozygous for a 
multicopy GFP transgene, as described previously [5]. All 
MommeD lines were maintained in this background and 
are homozygous for the GFP transgene. All experimental 
data: crosses to other lines, intercrosses, and so on, were 
carried out using heterozygous MommeD mice five gener- 
ations or more removed from the MommeD founder. 

Other mouse strains 

The congenic strain, LineSC, used for linkage studies, was 
produced by crossing LineS to C57BL/6J for 10 genera- 
tions, selecting for mice carrying the transgene by flow cy- 
tometry. Inbred C57BL/6J mice were purchased from 
ARC Perth (Perth, WA, Australia). Rifl^'^ mice were gen- 
erated by the Australian Phenomics Network from an em- 
bryonic stem cell clone carrying a trapped Rifl allele 
(A045A01; German Gene Trap Consortium), crossed to 
LineSC (C57BL/6J) and maintained on this background. 

Flow cytometry 

Mice were analyzed by flow cytometry at 3 weeks of age. 
A drop of blood was collected in Osmosol buffer and 
analyzed on a Guava easyCyte HT (Merck/Millipore, 
Darmstadt, Germany). The data were analyzed by using 
Guava InCyte software with a GFP-positive gate set to 
exclude 99.9% of wild-type erythrocytes. Histograms 
shown depict only the GFP fluorescence channel. 

Linkage analysis 

Heterozygous mutant MommeD mice, at least four 
generations down from the founder, were backcrossed 
twice to LineSC (see above) and phenotyped for GFP 
expression by flow cytometry. DNA from tail tips was 
used to perform linkage analysis. The Alumina GoldenGate 
genotyping assay (Mouse Medium Density Linkage Panel) 
was used for MommeD 13, MommeD 16 to MommeD 19, 
MommeD27, MommeDSO, MommeD32, MommeD34, 
MommeDSS, MommeD37, MommeD39, MommeD40 and 
MommeD42 on at least 12 wild-type and 12 heterozygous 
mice from each line. The Mouse Medium Density Linkage 



panel contains 766 measurable SNPs between C57BL/6J 
and FVB/NJ. Samples were genotyped following the 
Alumina protocol and genotype calls were made using the 
Genotyping module of the GenomeStudio vl.l software. 
Only samples with a call rate >95 were accepted. Linked in- 
tervals were identified based on LOD scores of 3 or higher 
(Additional file 7). 

For MommeDS, MommeD23, MommeD28, MommeD32, 
MommeD33 and MommeD36 the linked chromosomes 
were identified using microsatellite or SNP markers. Add- 
itional markers were used to reduce the linked intervals 
using wild-type and heterozygous mice. Mice wild-type for 
the mutation should only have C57BL/6J chromosomes in 
the linked region, while mice heterozygous for the mutation 
should carry both C57BL/6J and FVB/NJ chromosomes. 

For MommeD13, MommeDlG, MommeDlS, MommeD23, 
MommeD28, MommeD32, MommeD33 and MommeD35, 
mutations were found using Sanger sequencing. For the 
remaining MommeD lines deep sequencing approaches 
were used. 

Whole exome deep sequencing and Interval capture 

Mouse exomes were captured and sequenced using exome 
capture reagents from either Agilent (SureselectXT 
Mouse All Exon v.l, Santa Clara, CA, USA) or Roche 
NimbleGen (SeqCap EZ Mouse Exome, version Beta 2, 
110603_MM9_exome_rebal_2EZ_HXl, Madison, WI, 
USA). The Agilent captures were carried out as detaUed 
in the manufacturers protocol (version 1.1.1), with the 
foAowing exceptions: a Bioruptor (Diagenode, Liege, 
Belgium) was used for DNA shearing, producing frag- 
ment sizes of approximately 200 to 300 bp, and sample 
pooAng/multiplexing was omitted. The Bioruptor sonic- 
ation settings were 3 x 10 minutes on the low setting 
with cycles of 30 seconds on and 30 seconds off. The 
Roche NimbleGen captures were carried out as outlined 
in the Alumina optimized Roche NimbleGen SeqCap 
User s Guide (version 1.0), also using the Bioruptor for 
DNA fragmentation. 

For two mutant lines, MommeDS and MommeD34, a 
4.2 Mbp linked interval was sequenced in fuU using 
a custom capture array designed and manufactured 
by Roche NimbleGen. MommeDS was captured and 
sequenced by Roche NimbleGen using Roche 454 se- 
quencing (110317_MM9_mommeD8_Rocky_cap_HX3). 
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Three mice were sequenced; a wild type, a heterozygote 
and a homozygote. For MommeD34, a heterozygous mu- 
tant was captured and sequenced in-house using the 
same array capture design on 2.1M arrays with HX3 
mixers, according to the Roche NimbleGen Arrays 
Users Guide, lUumina Optimized protocol (version 1.0). 

SNP calling 

The sequencing reads from the targeted sequence cap- 
ture experiments were aligned to the mouse genome 
(build 37, mm9) using the program bwa version 0.6.1 
[57]. Datasets generated by Roche on the 454 platform 
were mapped using bwa bwasw [58] while the remaining 
datasets, generated on the Illumina HiSeq or GAIIx plat- 
forms, were mapped using bwa aln, with default settings, 
followed by bwa sampe with the default settings. The 
resulting sam files were converted to bam files and 
coordinate-sorted using SAMtools version 0.1.17 [59] 
and PGR duplicates were subsequently eliminated using 
the program Picard MarkDuplicates, version 1.48 [60]. 

For each sample, nucleotide variants were identified 
within the intervals for which linkage had previously been 
determined. This was achieved by creating a pileup file of 
the linked region using SAMtools mpileup, using the op- 
tion -q 20, followed by variant calling using the program 
Varscan version v2.2.8 [61] using the 'somatic' feature and 
the settings -min-coverage 15 and -min-var-freq 0.3. 
Varscan somatic calls sequence differences between a case 
and a control sample; as control the sequenced exome 
from a different ENU mutant was used. The case and con- 
trol exomes were constructed using the same library prep- 
aration methods and sequenced in the same deep 
sequencing run, but had different linked intervals. All 
exomes were mapped and processed in parallel, using 
identical settings, to minimize post-sequencing artifacts. 
The output from Varscan was manually screened for likely 
ENU mutations, appearing as heterozygous SNPs in the 
mutant and wild type in the control. These SNPs were in 
turn validated using the Sanger method. 

The custom capture and sequencing of the MommeD34 
linked region was carried out in-house without a matched 
wild-type control; instead, a merge of the three MommeDS 
deep sequencing samples previously sequenced by Roche 
was used as control. A merge of these 454 datasets was 
used in order to achieve greater read depth across the re- 
gion. SNP calling was subsequently carried out as de- 
scribed above. Exome sequencing datasets generated in 
this study are accessible via European Nucleotide Archive 
(ENA) under accession ERP003831. 

Genotyping 

Once the mutations had been identified, genotyping 
was carried out by either sequencing {MommeDlS, 
MommeD16, MommeDlV, MommeDlS, MommeD23, 



MommeD27, MommeD28, MommeDSO, MommeD32, 
MommeD34, MommeD35, MommeD36, MommeD37, 
MommeD39, MommeD40 and MommeD42) or, if the 
mutation had changed a restriction enzyme recognition 
site, by PGR and digestion and gel electrophoresis 
(MommeDS and MommeD19). Genotyping primers are 
provided in Additional file 8 and Sanger sequencing 
traces in Additional file 9. 

Embryo dissections 

All embryos were produced by natural matings and 
detection of a vaginal plug was counted as 0.5 dpc. 
Except where otherwise stated, embryos were produced 
by intercrosses. 

Introduction of the mutant lines into mice carrying the 
A'^y allele 

FVB/NJ mice heterozygous for the MommeDS mutation 
and homozygous for the GFP transgene were mated with 
G57BL/6J mice heterozygous for the A^^ allele. The coat 
color phenotype was classified at weaning by a trained 
observer as either yellow, mottled, or agouti. GFP ex- 
pression was determined by flow cytometry and used to 
classify the mice into mutants or wild-type for the 
MommeDS mutation. FVB/NJ mice carry the A locus 
and C57BL/6J mice carry the a locus. Yellow and mot- 
tled offspring carry the A^^ allele. All agouti-colored off- 
spring were genotyped by PGR, to assess whether they 
were A'^^/A and pseudoagouti, or A/a, as reported [56]. 

Bisulfite sequencing of the transgene HS-40 enhancer 

Bisulfite conversion of DNA was carried out using the 
EpiTect Bisulfite Kit (Qiagen, Doncaster, VIG, Australia) 
according to the manufacturer s instructions. At least two 
male adult spleens were used for each MommeD line. The 
bisulfite conversion rate was at least 97% and sequences 
were analyzed using the BiQ Analyser software [62]. Oli- 
gonucleotides to the bisulfite converted HS-40 enhancer 
region were as follows (5'-3'): GFPbisFl: AAAATAAAA 
TTTTTGGATTGTTATTATTATAA; GFPbisF2: ATATT 
TGTAATTTTAGTATTTTGGGAGGTT and GFPbisR: A 
ATGTGTAGTGAGTAGAAAGTGGATGTG. Gycling con- 
ditions were as follows: 94°G for 2 minutes for 1 cycle; 94°G 
for 30 seconds, 60°G for 30 seconds, 72°G for 45 seconds 
for 35 cycles and 72°G for 6 minutes for 1 cycle. 

Protein analysis 

Whole-cell extracts were prepared from various tissues of 
adult mice. Protein concentration was measured by BGA 
(Thermo Scientific, Waltham, MA, USA) and protein ly- 
sates were separated using polyacrylamide gels (Invitrogen, 
Garlsbad, GA, USA, BioRad, Hercules, GA, USA or EGL, 
Rydalmere, NSW, Australia). Antibodies used were as fol- 
lows: anti-y-tubulin (T5192, Sigma- Aldrich, St. Louis, MO, 
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X-axis represents the chromosomes and the y-axis is the LOD score. Peaks 
with a LOD score of 3 or higher are considered significant. 

Additional file 2: Flow cytometric profiles of the MommeDs. 

(a) Representative flow cytometry profiles show the percentage of GFP- 
expressing erythrocytes from wild-type and heterozygous mutant 
littermates (n = 3). The x-axis represents the erythrocyte fluorescence on 
a logarithmic scale and the y-axis is the number of cells detected at each 
fluorescence level, (b) Percentage of GFP-expressing cells in wild-type, 
heterozygous mutant and homozygous mutant (where viable) mice at 
three weeks of age (mean ± standard error of the mean). 

Additional file 3: Mapping interval for MommeDSO. (a) MommeD30 
was produced in the FVB/NJ strain of mice and mapped by crossing 
twice onto Line3C in a C57BL/6J background. The results of the 
genotype for seven SNP markers and one microsatellite marker 
surrounding the linked interval are shown. The number of mice classified 
into each haplotype is shown on top. Our estimate of the linked interval 
is between rs29539305 and rs33446195 on chromosome 17 (highlighted). 

(b) List of genes in the MommeD30 linked interval on chromosome 17. 
The Mus musculus EnsembI database (release 37) was used to export a 
list of transcripts (protein coding and non-coding RNAs) within the 1.9 
Mbp MommeDSO interval. 

Additional file 4: GFP expression in offspring of a Rifl^^ 
heterozygote crossed to LineSC A Rifl gene trap allele {Rifl^^) had a 
similar effect on transgene expression as that observed with the 
MommeDlS mutation, increasing the percentage of expressing cells in 
mice heterozygous for the gene-trap allele. 

Additional file 5: Novel mutant alleles of Bazlb, Smchdl, Trim28, 
Dnmtl and SmarcaS. (a) MommeD16 carries a point mutation resulting 
in a non-conservative amino acid change in the Wstf domain of Bazlb. 
(b) MommeD23 and MommeDSO carry point mutations in Smchdl. Both 
mutations introduce premature stop codons in the Smchdl protein, (c) 
MommeDSl carries a point mutation in Trim28 that results in an amino 
acid change in a highly conserved zinc finger domain, (d) MommeDS2 
carries a point mutation that results in an amino acid change in the BAH 
domain of Dnmtl. (e) MommeDSS and MommeDS? carry mutations that 
result in amino acid changes in highly conserved domains of the 
SmarcaS protein. 

Additional file 6: Embryonic development in MommeD mutants. 

Tabulated data shows the number of observed mice and in brackets the 
percentage of total, (a) Intercrosses of P6rm/^°"^"^^°^^ mice produce no 
homozygous offspring at three weeks, (b) srdl'^"^^^^^^ homozygous 
mice are embryonic lethal around 10.5 dpc. (c) Mice carrying mutations 
in 5^^39/7/""°"^"^"^'' are viable at three weeks, (d) uhrfl''"'''^'''^' 
homozygotes are embryonic lethal around 10.5 dpc. (e) Some 
homozygous Bozlb^"^"^^^^^ mice were obtained at three weeks but less 
than expected, (f) No {Smchdl''"'''''^''') or few {Smchdl'^"'''''^'') 
homozygous individuals were recovered from intercrosses. The 
Smchdl homozygotes that survived were males, (g) 
7r/m28^°"^"^^^"^' homozygous mice are embryonic lethal prior to 10.5 dpc. 
(h) Mice homozygous for the Dnmtl ^"'^'^''^^^ mutation die around 10.5 
dpc. (i) Intercrosses of Smarca5^°"^"^"°'^ and SmorcoS'^"'''"'^'' produced 
no homozygous offspring at three weeks. Timed matings show that 
survival of homozygous mutants at 12.5 to 14.5 dpc was rare in 
SmarcaS'^"'"'^'^'' mutants. In SmarcaS'^'''''^'^'' mutants homozygous 
embryonic death occurred prior to 14.5 dpc. 

Additional file 7: Linkage analysis using lllumina GoldenGate 
genotyping assay. 

Additional file 8: List of primer sequences. 
Additional file 9: Sanger traces. 



USA), anti-Rifl (abl3422, Abeam, Cambridge, ENG, UK), 
anti-Wiz (gift ft-om Yoichi Shinkai), anti-Smarccl (9053S, 
Cell Signaling, Danvers, MA, USA), anti-Rlf (abllSOll, 
Abeam - rabbit polyclonal) and (M05, clone2G2, Abnova - 
mouse monoclonal, Neihu District, Taipei City, Taiwan), 
anti-Pbrml (Rabbit polyclonal ABE70, Millipore), anti- 
Suv39hl (clone D11B5, Cell Signaling), anti-Uhrfl 
(abl51187 Abeam), anti-GAPDH (D16H11 XP, Cell Signal- 
ing) and anti-Snf2H (ab3749. Abeam). 

Cell fractionation 

HeLa cells were trypsinized, pelleted and washed twice with 
phosphate-buffered saline before re-suspending in ice-cold 
cytoplasmic extraction buffer (10 mM Tris pH 7.4, 10 mM 
NaCl, 3 mM MgCb, 0.1% NP-40). A small aliquot of this 
suspension, the total cell lysate, was transferred to a sepa- 
rated tube. The remaining suspension was then centriftiged 
(500 rpm, 10 minutes, 4°C) to pellet the nuclei. The super- 
natant (cytoplasmic fraction) was aliquoted into a separate 
tube. Nuclei were re-suspended in an 8 M urea lysis buffer 
(8 M urea, 1/10 vol. glycerol, 1/20 vol. 20% SDS, 1/2,000 
vol. 1 M dithiothreitol, 1/100 vol. 1 M Tris, pH 6.8). 

RNA isolation, cDNA analysis and quantitative real-time 
RT-PCR 

Total RNA was extracted from various tissues using TRI re- 
agent (Invitrogen). cDNA was synthesized from total RNA 
using SuperScriptlll reverse transcriptase (Invitrogen) or 
AMV reverse transcriptase (Roche) and random hexamer 
primers. Quantitative real-time PCR was performed with 
the Platinum SYBR Green qPCR Super Mix -UDG 
(Invitrogen) with primers designed to span exon/intron 
boundaries. All reactions were performed in triplicates and 
normalized to Hprt or Gapdh. PCRs were run on a Viia7 
(Applied Biosystems, Mulgrave, VIC, Australia) or on a 
Corbett Research Rotor-Gene (Qiagen). Cryp-Skip [63] was 
used for splice site prediction. cDNA from mutant alleles 
was sequenced using Sanger sequencing. Primer sequences 
are provided in Additional file 8. 

Statistical analysis 

Statistical significance of quantitative data was deter- 
mined by two-tailed Students T-test. F-test was used to 
test whether variance was significantly different between 
wild-type and mutant groups. The proportions of geno- 
types were compared to expected Mendelian ratios using 
a test. For all datasets a minimum of three biological 
replicates were analyzed. 

Additional files 



Additional file 1: Linked intervals. Manhattan plots showing linked 
intervals identified by lllumina GoldenGate SNP genotyping analysis. The 



Abbreviations 

bp: base pair; dpc: days post-coitum; ENU: N-ethyl-N-nitrosourea; 
FSHD2: facioscapulohumeral dystrophy type 2; GFP: green fluorescent 
protein; H3K9: histone H3 lysine 9; Momme: Modifier of murine metastable 
epiallele; PCR: polymerase chain reaction; SNP: single-nucleotide 
polymorphism. 
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